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Abstract 

The extent to which information about fluctuations in hadron-nucleon total 
cross sections in the frozen approximation can be extracted from very high 
energy hadron-nucleus total cross section measurements for a range of heavy 
nuclei is discussed. The corrections to the predictions of Glauber theory due 
to these fluctuations are calculated for several models for the distribution 
functions, and differences of the order of 50 mb are found for heavy nuclei. 
The generating function for the moments of the hadron-nucleon cross section 
distributions can be approximately determined from the derivatives of the 
hadron-nucleus total cross sections with respect to the nuclear geometric cross 
section. The argument of the generating function, however, it limited to the 
maximum value of a dimensionless thickness function obtained at zero impact 
parameter for the heaviest nuclear targets: about 1.8 for pions and 3.0 for 
nucleons. 
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I. INTRODUCTION 



There has recently been a revival of interest in the use of the cross section distribution 
function to describe certain features of hadronic collisions |l]-[3] . The general approach, which 
goes back to the work of Good and Walker [ffl in 1960, begins with the composite nature 
of hadrons and assumes that at high energies the internal degrees of freedom are "frozen" 
during the collision M. Many features of high energy reactions can then be discussed, at 
least approximately, in terms of a single distribution function which gives the probability 
that the initial hadron is in one of the configurations which interacts with a given total cross 
section. Discussions of this distribution function for incident pions and nucleons interacting 
with nucleons, including estimates of the 2nd and 3rd moments, are given in Refs. [1] and 
[2]. The discussion below tries to answer the question of whether more information about 
these distributions, and thus more information about the composite nature of the hadrons, 
can be obtained by a careful study of hadron-nucleus total cross sections. 

The answer to this question is closely related to the importance of inelastic intermediate 
states in hadron-nucleus scattering, which has been discussed extensively 0. In particular 
the systematic experimental study of Murthy et al. || confirmed earlier theoretical estimates 
of these effects, including their energy dependence. Here a similar discussion is presented in 
the language of the cross section distribution function. This approach is valid only in the 
high energy limit, where the frozen approximation has a chance of being valid, but has the 
advantage of including all orders of inelastic scattering, which might well be important for 
heavy nuclei. 

Section II gives a brief review of how the distribution function is defined, and what is 
known about it. Then Sec. Ill summarizes the approximations and assumptions which are 
necessary to obtain a simple formula for the total hadron-nucleus cross section in terms of 
the generating function for the reduced distribution function, the argument being a reduced 
nuclear thickness function. In Sec. IV the limiting case of uniform nuclear density is dis- 
cussed. It is shown that in this case the generating function of the distribution function is 
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just the derivative of the total cross section with respect to twice the cross sectional area 
of the nucleus. This result shows clearly that the amount of information we can hope to 
obtain from nuclear cross sections is limited because the argument of the generating function 
is proportional to the nuclear thickness function, and thus limited by the limited sizes of 
stable nuclei. This discussion is extended to the more realistic case of Woods-Saxon nu- 
clear densities in Sec. V. For these densities the derivatives mentioned seem to approach 
the generating functions from above for large nuclei. In addition the numerical results of 
this section give some idea of the sensitivity of these quantities to reasonable changes in 
the distribution function, and thus an estimate of the experimental accuracy required to 
distinguish among different possibilities. These results are discussed in Sec. VI, along with 
suggestions for further investigations in this area. 



II. CROSS SECTION DISTRIBUTION FUNCTION 

At high energy it is assumed that the internal configuration, here labeled by a Greek 
letter a, (3, etc., of a hadron interacting with a target is "frozen" during the interaction so 
that the transition amplitudes are diagonal in these states: 

<J3 | F op | a) = 5p a F a , (1) 

where F op is the transition operator in the space of the internal degrees of freedom of the 
incident hadron. The amplitudes for transitions between ordinary hadronic mass eigenstates 
i and j on a given target are then 

(j \Fop\i) = £ (j | a)F a (a | i), (2) 

a 

where i and j must represent states with the same flavor quantum numbers, so j must be 
either i or a state which can be diffractively excited from it. Using the same normalization 
as Blattel et al. |],|2| the total cross section for hadron i incident on the target is 

Oi = 47clm(i | Fop | i) = £ a a \ (a \ i) \ 2 , (3) 

a 
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where a a is the total cross section for the hadron in configuration a to interact with the 
target. 

If the cross section distribution function Pj is defined as 

Pi(a) = £ I H *> \ 2 6(a-a a ), (4) 

a 

this becomes 

c?i = (« | o" op | *) = / -Pi(cr) cr da, (5) 



i.e. (7j is just the average of a over the distribution Pj. The same distribution can also be 
used to calculate the forward total differential cross section for diffractive scattering, if the 
transition amplitudes are assumed to be pure imaginary, for then 



dt 



£ I 0' I F °p I *) 

_2\ /_\2 



«Oi " (*>?) / (167T). (6) 



For scattering in the high energy regime discussed by Blattel et al. |T|fJ], (cr) ~ 24mb for 
pions on nucleons while (cr) ~ 40mb for nucleons on nucleons. Experimental results on high 
energy forward diffraction dissociation on nucleons gives (cr 2 ) m 1.25 (cr) 2 for nucleons, while 
(cr 2 ) m (1.4 - 1.5) (cr) 2 for pions. It is also possible to put constraints on the 3rd moments 
of the distributions using diffraction dissociation on deuterons, but only by making fairly 
strong assumptions about the transitions between excited states. 

For the developments below it is convenient to introduce reduced distribution functions 

/,(*) ee (a) P t ((a) x), (7) 

so that 

(x n )i ee J dx h{x) x n = (a n ) t I (cr)". (8) 
The simple Poisson-like distributions 

f n (x) = [(n + l) n+1 / n\}x n exp [-(n + l)x] (9) 



will be used below for illustrative purposes because they lead to analytic formulas in some 
cases. The generating function, for example, is 

71+1 

(10) 

leading to the moments 

(o. = ■ (id 

n\ (n + m) m 

Blattel et al. suggest several different forms for the pion-nucleon and nucleon-nucleon 
cross section distribution functions, all of which fall off with x more rapidly than the / n 's. 
Based on the expected behavior of the pion and nucleon internal wave functions they also 
assume that f(0) vanishes for nucleons, but not for pions . The behavior at large x, however, 
is more problematic. The main object of this paper is to see if a careful and systematic study 
of hadron-nucleus total cross section can shed any new light on these distribution functions. 



(exp {-xt)) n 



(n + 1) 



(n + l + t) 



III. HADRON-NUCLEUS CROSS SECTIONS 

In Glauber theory, which itself uses the frozen approximation for the nucleons in the 
nucleus, the total cross section for scattering of a hadron % from a nucleus A can be written 

as 

<n(A) = 2 J d%Y i {A ) b) 1 (12) 

where 

ri(A,b) = l-(A\U a (1 -7<,«(6)) I A), (13) 

with 7i jQ , the analogue of T for the scattering of the hadron i from the ath nucleon in 
the nucleus. In "standard" Glauber theory, where the internal degrees of freedom of the 
hadron are ignored, 7; iQ is just a numerical function of b and r Q , the position of the ath 
nucleon (except for possible spin- and isospin-dependence, which we ignore here). In general, 
however, 7 should be interpreted as an operator in the internal space of the hadron, and 
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Ti(A, b) = 1 - (A | {i | U a (1 - 7aop (6)) | i) | A). (14) 

This expression can be simplified considerably if a number of assumptions are made. 
First, if the nuclear wave function is completely uncorrelated 

T t (A,b) = 1 - (i | [1 - (A | 7lap (6) | | ^). (15) 

Then, if the range of 7 is much less than the nuclear size, 

T t (A,b) = l-(i\[l-^T(b)] A \i), (16) 

where T(b), the thickness function, is the integral of the nuclear density along a straight line 
at impact parameter b, and <j op is twice the integral of 7op over impact parameter space. In 
general 7 op and a op are both operators in the internal space of the hadron taking complex 
values, but if the imaginary part is ignored then <r op is exactly the operator described by the 
distribution function Pj introduced in Sec. II. Finally, if A is large and <T op T(b) takes only 
small values, then 

r,(A,b) = l-(i|exp[-^T(b)A]|<> 

= 1- (i I exp [-t (A,b)x op ] I (17) 

where 

t(A,b) = AT (A, b)(a) / 2 (18) 
is the dimensionless reduced thickness function and 

Xop = Pop I (c) (19) 

is the dimensionless reduced cross section operator. In other words, is determined by 

(i I exp (—t x op ) I i) — J dxfi(x) exp (—t x), (20) 
the generating function for the reduced distribution function /j(x). 
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The reduced thickness function t(A,b) is expected to be largest for b=0 and to increase 
nearly monotonically with A. For the largest stable nuclei (A=238, say) , t max 3.0 for 
incident nucleons and 1.8 for pions. We can therefore not hope to get information on the 
generating function for t greater than these values from nuclear total cross sections. 

If fi(x) does not converge rapidly enough at large x then the short-range and exponential 
approximations may be invalid. The short range approximation is essential if the nuclear 
cross sections are to be expressed in terms of /i(x), but the exponential approximation is 
to some extent just a convenience: a version of most of the results below could be obtained 
without it. For small t 

(i | exp (-t x op ) | i) = 1 - t + t 2 (i | x 2 | i) I 2 - ... , (21) 

provided the series converges. (It does converge for all t for all the examples of Refs. [1] and 
[2], but for the Poisson-like distributions the radius of convergence in n+1.) Since (i | x 2 | i) 
is about 1.4 to 1.5 for pions and 1.25 for nucleons, and some constraints on (i | x 3 | i) can 
be obtained from forward diffractive cross sections on deuterons, the generating function is 
already fairly well determined for small t. Nuclear cross sections will provide new information 
only if they can be used to constrain the generating function for values of t larger than about 
1. 

The behavior of the generating function for very large t is determined by the behavior 
of the distribution function near x=0: 

(i | exp (-t x op ) | i) — + - + ... , (22) 

where f^(x) is the mth derivative of fi(x). If the generating function could be determined 
for large enough t, then, one could check the assumptions of Refs. [1] and [2] that fi(0) 
vanishes for nucleons, but not for pions. Unfortunately, as noted above, the range of t is 
limited by the limit on the sizes of stable nuclei to values less than 3. 

For many purposes it is useful to decompose Eq. (12), together with Eq. (17), as 

a,(A) = af\A)-of\A), (23) 
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where 

af\A) = 2 J d 2 b G {G) (t(A,b)) (24) 

is the Glauber result, and 

a\ D \A) = 2 J d 2 bG {D) (t(A,b)). (25) 

is the decrease due to the dispersion in the hadron-nucleon total cross section. (It is easy to 
show that G(°) cannot be negative and vanishes only if there are no fluctuations in a.) The 
integrand in Eq. (24) is simply 

(t) = 1 - exp (-*), (26) 

while 

G {D) (t) = (i | exp (-t x op ) | i) - exp (-*). (27) 

Because is a fairly small correction to <7i(A), in applying these formulas to ex- 
perimental data it may be sufficient to include non- dispersive corrections, such a nuclear 
correlations and corrections to the short range approximation, only in <j( G \ leaving cr(D) as 
the simple expression above. 

IV. UNIFORM DENSITY LIMIT 

As noted above, new information about the generating function can be obtained from 
total cross sections only for heavy nuclei. The nuclear densities p(r) are then nearly constant 
in the nuclear interior, falling quite rapidly to zero at the nuclear surface. In this section 
the extreme uniform density limit is assumed: 

e (r) = 6(R-r) I (AtvR 3 / 3), (28) 

where R ~ r A 1//3 , with r ~ 1.1 fm, is the nuclear radius. In this limit, replacing A by the 
equivalent and more convenient label R, 
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t(R,b) = t(R,0)yfl - (b/R) 2 9 (R-b), (29) 

where 



t(R,0) = aR, 



(30) 



with 



a = 



a I {Aixrl I 3). 



(31) 



With this simple expression for t(R,b), for any expression for a contribution to a total 
cross section of the form 



In other words, the slope of the cross section contribution as a function of twice the 
nuclear geometric cross section is just the function G evaluated at the maximum value of 
t(R,b) for that nucleus, obtained at zero impact parameter. Since G can be chosen to be 
1 minus the generating function, or the generating function minus exp(-t), it is then trivial 
to extract the generating function at t(R,0) from the derivative of the cross section. This 
relation, however, is exact only for the uniform density case. In the next section the accuracy 
of this relation for more realistic Woods-Saxon densities is studied. 



A more realistic nuclear density, especially for the heavy nuclei of main concern here, 
has the Woods-Saxon form 




(32) 



it is easy to show that the derivative 




(33) 



V. WOODS-SAXON DENSITIES 



g(R,r) 



Qo (R) 



(34) 



1 + exp ((r - R) / a ) ' 
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where ao ~0.523 fm is the surface thickness parameter and 



1 / g (R) 




2 _^0_ 



(35) 



is chosen so that the volume integral of p is unity. (The parameter R for a number of different 
nuclei have been determined by fits to total cross sections at energies below 30 GeV/c, where 
inelastic contributions are expected to be small ||[7|- More exact nuclear densities can be 
determined from electromagnetic form factors, but the Woods-Saxon form is sufficient for 
present purposes.) Although there is an exact analytic expression for 



the thickness function for other values of b must be obtained by numerical integration. For 
this density the derivative relationship found in Sec. IV is not exact, although one might 
expect it to be a good approximation for large R where the density is quite uniform over 
most of the nucleus. In this section numerical results determining its accuracy are presented. 

Fig. 1 shows er D (R) as a function of 27rR 2 for 4 different pion-nucleon cross section 
distribution functions for pion- nucleus scattering. Distributions (b), (c), and (d) are the 
same as in Fig. 1 of Ref . [2] , while distribution (a) is the superposition of two Poisson-like 
distributions introduced in Sec. II. Fig. 2 shows the analogous curves for nucleon-nucleus 
scattering, together with two experimental points, along with their error bars, from Ref. [6] . 
Here curve (a) is a two-term exponential while (b) and (c) are the n=2 and n=6 distributions 
shown in Fig. 4 of Ref. [1]. (The n=10 distribution gives results which are almost identical 
to those for n=6. In using the distributions from Refs. [1] and [2] we have adjusted the 
parameters slightly so that (1) = (x) = 1.0 more precisely.) For small 27rR 2 all the 
distributions give nearly the same cx D (R) because only low moments contribute, but for 
large 27rR 2 there are differences of order 50 mb. Since these are of the same order as the 
experimental errors in Ref. [6] , it is clear that more accurate experiments on heavy nuclei at 
very high energy are needed if the hadron-nucleus cross sections are to discriminate among 
different cross section distribution functions. 



t(R, 0) 



Aag {R + In [1 + exp (-R / a ))} , 



(36) 
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In Fig. 3 the slopes in Fig. 1 are compared with the G D (t(R,0))'s for the various distri- 
butions. At large R the two quantities are similar, but G D is uniformly below the derivative, 
although the difference slowly decreases with increasing R. The corresponding results for 
nucleon-nucleus scattering are shown in Fig. 4. It seems that one might get a reasonable 
estimate for G D if o" D (R) could be measured accurately enough. This would be particularly 
true if an estimate for the surface corrections to the derivative relationship at large R could 
be found. 

VI. CONCLUSION 

The results above suggest that accurate measurements of the total cross sections for the 
interactions of hadrons with a range of heavy nuclei at very high energies could increase 
our knowledge of the hadron-nucleon total cross section distribution functions, fixing the 
generating functions up to fairly large arguments and selecting among models suggested 
by theoretical prejudices and constraints on the first few moments. Very high energies 
are required for the validity of the "frozen configuration" approximation, while errors in 
measuring cross sections ranging to more than 3000 mb must be restricted to perhaps less 
than 10 mb. The dispersive effects are larger for pions than for nucleons, but unfortunately 
the maximum value of the argument of the generating function is smaller for pions because 
of their smaller cross sections on nucleons. 

Even if the results of such difficult measurements were available, however, more work 
is required before they could be unambiguously interpreted in terms of the distribution 
functions. The formulas above depend upon a number of simplifying approximations: the 
frozen approximation, the neglect of nuclear correlations, the short range approximation, and 
the neglect of the imaginary part of the hadron-nucleon scattering amplitude. Fortunately, 
as noted in Ref. [6], the energy dependence of the total cross sections gives us a handle on 
many of these. In the range of energies high enough that the Glauber approximation is 
valid but low enough that inelastic intermediate states do not yet contribute significantly, 
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the nuclear total cross sections (as well as the differential cross sections) can be used to 
check assumptions about the nuclear wave functions. At the other extreme, the total cross 
sections at very high energies should approach constant values (aside from slow changes due 
to the energy dependence of the hadron-nucleon parameters) if the frozen approximation is 
valid. 

It would be useful to have more information about the surface corrections to the relation 
between the generating function and the derivative of the total cross section as a function 
of R derived in Sec. IV in the uniform density approximation. Preliminary investigations 
suggest that for Woods-Saxon distributions these corrections should decrease as an inverse 
power of R. If this is true it should be possible to obtain more accurate values for the 
generating functions by studying the dependence of the total cross sections on R. 
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FIGURES 

FIG. 1. Dispersive contribution to pion-nucleus scattering as a function of 27rR 2 , where R is the 
radius parameter in the Woods-Saxon nuclear density. Curves b, c, and d are calculated using the 
corresponding distribution functions shown in Fig. 1 of Ref . [2] , while curve a is calculated using a 
distribution function which is a superposition of the n=0 and n=l exponential distributions defined 
in Sec. II. This distribution was chosen to match distribution c at a=0 and has (a 2 ) / (a) 2 =1.7 
instead of the value 1.5 of distributions b,c,and d. 

FIG. 2. Dispersive contribution to nucleon-nucleus scattering as a function of 27rR 2 , where R is 
the radius parameter in the Woods-Saxon nuclear density. Curves b, and c are calculated using the 
n=2 and n=6 distribution functions shown in Fig. 4 of Ref. [1], while curve a is calculated using a 
distribution function which is a superposition of n=l and n=4 exponential distributions chosen to 
also give (a 2 ) = 1.25 (cr) 2 . All these distributions vanish at a=0. The experimental points with 
their errors are taken from Ref. [6]. 

FIG. 3. The solid curves are the derivatives of the corresponding curves in Fig. 1, while the 
adjacent dashed curves are the G( D ) of Eq. (27) for the same distributions, evaluated at i(R,0) as 
given by Eq. (30). 

FIG. 4. The solid curves are the derivatives of the corresponding curves in Fig. 2, while the 
adjacent dashed curves are the G( D ) of Eq. (27) for the same distributions, evaluated at i(R,0) as 
given by Eq. (30). 
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